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Abstract 

Background: A successful development of herbivorous insects into plant tissues depends on coordination of 
metabolic processes. Plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense 
response. To understand the transcriptional basis of this response, we compare gene expression profiles of two coffee 
genotypes, susceptible and resistant to leaf miner {Leucoptero coffello). A total of 22000 EST sequences from the Coffee 
Genome Database were selected for a microarray analysis. Fluorescence probes were synthesized using mRNA from 
the infested and non-infested coffee plants. Array hybridization, scanning and data normalization were performed 
using Nimble Scan® e ArrayStar® platforms. Genes with foldchange values +/-2 were considered differentially expressed. 
A validation of 18 differentially expressed genes was performed in infected plants using qRT-PCR approach. 

Results: The microarray analysis indicated that resistant plants differ in gene expression profile. We identified relevant 
transcriptional changes in defense strategies before insect attack. Expression changes (>2.00-fold) were found in 
resistant plants for 2137 genes (1266 up-regulated and 873 down-regulated). Up-regulated genes include those 
responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and 
maintenance. Also, our analyses indicated that differential expression profiles between resistant and susceptible 
genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant 
plants, as a priming mechanism. Validation of selected genes pointed to four selected genes as suitable candidates 
for markers in assisted-selection of novel cultivars. 

Conclusions: Our results show evidences that coffee defense responses against leaf-miner attack are balanced 
with other cellular functions. Also analyses suggest a major metabolic reconfiguration that highlights the complexity 
of this response. 
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Background 

Once a plant recognizes a pathogen attack, the metabolism 
must balance demands for resources to support defense 
versus requirements for cellular maintenance, growth and 
reproduction [1-4], Defense mechanisms involve a shift on 
metabolism, activating specific pathways such as synthesis 
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of secondary metabolites, programmed cell death, ions 
translocation. Concomitantly, can occur a shutdown of 
other metabolic pathways not directly involved with 
defense response, such as those associated with growth 
and reproduction, A resistance and/or tolerance trait is 
attributed whenever this defense response is successful 
in controlling pathogen or herbivore attack. 

Genetic control of metabolic re-programming is normally 
triggered by few resistance genes which are seek out to be 
transferred to other plants. However, resistance response 
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involves changes in other genes, not usually identified, 
with determinant roles in the overall response. Therefore 
for an effective transference of resistance traits, to know 
how these genes interact during re-programming of plant 
metabolism is essential 

Among available methods for high-throughput analysis 
the microarray is a powerful tool for large-scale gene 
expression studies in many plant species with whole 
genome sequenced: potato [5,6], tomato [7,8], soybean 
[9,10], wheat [11], barley [12,13] maize [14,15], grape [16], 
pine [17], Arabidopsis [18-20]. The main advantage of 
microarray analyses is to evaluate the expression of large 
number of genes in different genotypes, organs, tissues, 
treatments, using the same set of genes. These genes 
can be compared during different biological situations 
allowing both an association with metabolic pathways and 
establishement of their role on resistance response. Several 
studies have been carried out using microarray analysis 
to identify genes associated with plant defense [21-25]. 

The leaf-miner, Leucoptem coffeella (Guerin-Meneville, 
1842) (Lepidoptera-Lyonetiidae) is a specialist parasite 
of Coffea species. Upon oviposition on leaves, ecloded 
larvae feed directly from parenchyma tissues, leading to 
a reduction of foliar surface and eventual leaf drop 
[26]. This damage results in reduction of photosynthetic 
area and plant survival. In Brazilian breeding programs 
resistance genes from C. mcemosa have been transferred 
to the susceptible C. arabica by controlled crosses, and so 
far a large number of hybrid progenies are under selection 
for resistance to leaf-miner [27]. Although defense mech- 
anism to leaf-miner is not understood yet, genetic analysis 
demonstrated that resistance to the insect is dominant 
and controlled by two complimentary genes [28]. At 
the molecular level, there is little information regarding 
gene expression on coffee plants during defense response. 
Using subtractive hybridization methodology (SSH), 
Mondego et al. [29] found differentially expressed genes 
in coffee plants upon leaf-miner infestation, among which 
a miraculin-like encoding gene was significantly overex- 
pressed in resistant coffee plants. Differential expression 
of defense-related genes such as lipoxygenase, glutathione 
transferase, protein-kinase receptor and glucanase was 
observed in response to leaf-miner infestation [30]. 
However, the expression profiles indicate that differences 
results from gene expression timing along insect infection 
rather than with gene regulation. 

Despite the efforts of breeding programs to develop novel 
coffee cultivars bearing leaf-miner resistance, selection 
of progenies homozygotes for this trait is difficult [27], 
as advanced generations are still producing susceptible 
plants. Therefore, information regarding molecular control 
of resistance response as well as identification of candidate- 
genes associated with these processes will contribute with 
assisted-selection. 



In this context, the aim of our study was to explore 
transcriptomic differences throughout insect infestation, 
in susceptible and resistant C. arabica plants challenged 
by L, coffeella, using microarray technology. The arrays 
were developed using coffee -specific oligoprobes designed 
based on gene sequences available at the Brazilian Coffee 
Genome Project [31]. The database contains a collection 
of around 32,000 gene sequences, covering most of the 
C. arabica genome [32]. Besides this, we selected a group 
of candidate-genes to be used as molecular markers for 
assisted-selection. As far as we know, this is the first 
report of a large-scale transcriptional profile analysis 
used to study gene expression changes in coffee plants 
in the presence of an herbivore insect. 

Results 

Microarray analyses were performed to characterize large- 
scale gene expression profiles during leaf-miner develop- 
ment on coffee leaves. The analyses included a hybridization 
of a 135 K array with 6 different samples, corresponding 
to time-course infestation stages in both resistant and sus- 
ceptible plants. The arrays contain sequences of around 
33 K genes identified in EST libraries prepared from 
different physiological and metabolic situations [31]. A 
minimum of 6 - 24mer match probes for each selected 
gene were used for the array set up. The arrays were 
hybridized with probes corresponding to 3 treatments of 
both susceptible (S) and resistant (R) leaves: non-infestated 
(TO), after oviposition and egg-eclosion (Tl) and damaged 
by insect feeding (T2). 

Initially, differential expression patterns were identified 
using statistical analysis, and specific transcriptional profiles 
were established for each evaluated interaction. In a second 
approach, genes exhibiting differential expression among 
genotypes and treatments were submitted to in silico eval- 
uations to classify and categorize those genes regarding 
their possible molecular functions and metabolic pathways. 
Finally, a group of 19 genes involved with defense-related 
mechanisms, exhibiting regulated expression, were further 
characterized using qRT-PCR. 

Microarray and statistical analysis 

A total of 2141, 2359 and 2257 differentially expressed 
genes were identified from in silico analyses of raw hybrid- 
ization data considering 3 interactions: TOR X TOS, TIR X 
TIS and T2R X T2S (Table 1). Comparing TOR X TOS 
we observed higher differential expression levels where 
foldchange values range from 400 and 1000 times in 
up-regulated genes, and from 200 to 400 times in 
down-regulated genes in resistant leaves. The other 
interactions exhibited foldchange values ranging from 
150 and 350 times in down-regulated, and from 10 to 
15 times in up-regulated genes (Table 2). 
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Table 1 Distribution of regulated genes in each 



interaction analyzed 



Genes 


TOR_TOS 


T1R_T1S 


T2R_T2S 




(non infested 
plants) 


(Egg hatching) 


(Egg edosion) 


Up regulated 


1268 


1231 


889 


Down regulated 


873 


1128 


1368 


No differences 


19057 


18837 


18939 



Different regulation profiles for defense response were 
observed among differentially expressed genes. The first 
group includes chitinase and polygalacturonase genes, 
regulated differently in resistant and susceptible leaves 
from TO to T2. The second group includes genes that were 
up or down regulated in control resistant or susceptible 
leaves, but throughout insect development inverted their 
expression pattern. Example of this group is the gene 
encoding the enzyme polyphenol oxidase, up regulated 
at TO in resistant plants and down regulated at Tl and T2. 
And a third group includes genes that were up and/or 
down regulated in response to the leaf-miner infestation, 
either in susceptible or resistant leaves. 

Interaction between resistant and susceptible genotypes 
without insect infestation (TOR x TOS) 

Since higher values for differential gene expression were 
observed when comparing TOR X TOS samples, we chose 
this interaction for further analyses and selection of 
candidate-genes for validation. In this interaction, 2141 
genes exhibited differential expression, 1268 were up 
regulated and 873 were down regulated. 

Regulated genes from the TO samples were functionally 
characterized into three gene ontology categories - 
molecular function, component cellular and biological 
function - and grouped according their metabolic categor- 
ies (Figures 1 and 2). A description of the first one hundred 
most variable genes, both up and down regulated, is shown 
on Tables 3 and 4 and illustrated in Figures 3 and 4. Several 



Table 2 Survey of differentially expressed genes in 
all considered interactions including number and 
fold-change values range 



TOR_TOS 


T1R_T1S 


T2R_T2S 


Range of "fold change' 


14 


9 


10 


1000-100 


32 


17 


23 


99-40 


57 


51 


39 


39-20 


103 


96 


125 


19-10 


1062 


1058 


1171 


9-2 


858 


1117 


880 


(-)2- (-)9 


15 


11 


9 


(-)IO- (-)499 


Total 2141 


Total 2359 


Total 2257 





contigs had no correspondence with defined categories. 
A larger number of genes are up-regulated rather than 
down-regulated, and associated with primary and cellular 
metabolism, which included functions such as ion, protein 
and nucleic acid binding, hydrolase and transferase activ- 
ities, among others. As these differences were observed 
in non-infested leaves, possibly a different transcriptional 
programming takes place in resistant leaves. This may 
result in a pre-defense status, enabling resistant plants 
to a faster defense response upon leaf-miner attack. 

Metabolic pathways 

The categorization of annotated genes using Kegg database 
was performed with the first 100 up and down regulated 
genes (Figures 1 and 2). Most of these genes are from 
starch and purine metabolism, and several others are 
involved in primary metabolism. Three main metabolic 
pathways are highlighted here: citrate metabolism, linoleic 
acid metabolism and phenylpropanoids metabolism. 

We choose the citric acid cycle for further character- 
ization as previous analyses using NMR indicated that 
lower levels of malate, a metabolite resulting from the 
citrate metabolism, are observed in resistant coffee leaves 
[33]. Several genes encoding citric acid cycle enzymes 
exhibit differential-expression (Figure 5). Expression of 
isocitrate lyase gene is repressed in resistant genotypes 
at TO (fold change value -8,84), suggesting that synthesis 
of malate may be deficient, and therefore low levels of 
malate may accumulate in those leaves. However, this gene 
is up-regulated upon oviposition and egg ecclosion (fold 
change value 2,33). 

The linoleic acid is the first substrate of the Jasmonic 
acid (JA) pathway, a major signaling pathway during 
herbivore-defense responses. Control resistant plants 
(TO) show up-regulation of jasmonate 0-methyltransferase 
and lipoxygenase while differential expression for these 
genes was not observed at any stage in susceptible 
genotypes (Figure 6). Also, 13 genes from the alpha- 
linoleic acid metabolism and 57 genes from jasmonate 
biosynthesis were regulated in resistant plants (Figure 6). 
For instance, transcripts of enoyl-CoA hydratase and 
phospholipase A2 were four times more expressed in 
TO than in Tl in resistant genotypes (Figure 6), but 
increased only at later stages in susceptible leaves. 

Phenylpropanoids are major plant phtytoalexins, part of 
the secondary metabolism (Figure 7). Twenty-seven genes 
from phenylpropanoids synthesis exhibited differential 
expression at TO, with foldchange values ranging from 
9 to -5. Transcript levels of phenylalanine ammonia 
lyase (PAL), the enzyme that catalyzes the first step of 
the pathway, is up-regulated only at TO in resistant plants 
(2,05), and this level decreases along insect development. 
In susceptible plants, PAL transcript levels increase at final 
steps of insect infection, T2 (2,55). However, genes 
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KEGG pathways (up-regutated genes) 



Stireh and sucios« mslabolism 
Pyrirts mslabo^ism 
Amtno sugar and nuclMlide sugar metabolism 
iMethan^ melaboli^iTi 
alpha-Linolemc acid melabolism 
NiUcfjdrii melabolisin 
Glyoxylale and dicarboxylata m^Eabolism 
AscDrbat« and aldaratfr matabohsin 
Amirtdbcin££]at& d&gradalicin 
FlavonO'id biosynthe'Sls 
Glycolysis / GiuconeogenMift 
Tryptophan rvielabolisfin 
Mfllabolisin fsf x«nobiolics by cylo&hfom^ P450 
Ai^imrt^ and piolinB m^taboli^m 
Ponlose and ^lucuronale intfircortv^r^i^n^ 
T cell (Qceptar signaling patrtway 
Pn^nylalamrt^ melabolisitJ 
b«ta-Alanir»a melabotism 
Glylglfiione melabolisrn 
Pentose phosphate pathway 
Oxidative phDsphorylalian 
Oiug melabcli^m - GyEochiame Pi5<l 
Linclfiic acid melabcli^iiin 
Galactose meiaboli^in 
PnppanoatB melaboii^m 
Riboflavin! melaboli$ni 
Cyanoammo acid nrtalaboli^m 
Glycofolipid nrtelaboli«m 
Fatty acid nislabo4ism 
C«[tH>n rtxation in pholosynlh^ic oi^ani^ms 
eiosynttie^is of ijn&atuial;ed TaCly acids 
Tyrosine; melabdism 
Phenyl pi^panoid bipsynlbesis 
QutanoatB nielaboHism 
Qiycinv, sanii^ and ihiBonme m&tsboiism 
Glycenjpliospbolipid meiaboJism 
Pofphyiin and chloiqpbyll m diabolism 
TrOppne, pi pyridine and pyridine; alkaloid bioEynlhesis, 
Lysine degradation 
Tcrponoid backbone biDsyntbesi^ 
Limon^n^ and pmeng' degradation 
V*lin«. laucm^and isoteitiicin^ d^^radaliion 
StreploiTiycin biosynibesis 
Pynjvete meiabclism 
Fnjctos^ and mannose metabQlism 
Sen^oate degradatioHi 
□ileipenoid biosynlbe$i$ 
Gemniol degradatiwi 
□lug metabolism - other enzymes 
Pyninidinie metabolisni 
Capral^cta m degradatiofi 
Plavane and llavonal biDsynlhesis 
Fatty acid elongation 
Vitamin 86 melaboHism 
l$oquinolina alkaloid biosynlbasifi 
Inosliot phosipriate meiabdism 
Cystaino and melbionino melabolism 
Suff Lir melabolism 
Alanlna, aspartat* and glutamata malabolism 
Anihocyanih biosynltio$i$ 
Fatty acid biosyntha«i« 
Phanylalanlne. lyrosme and tryptophan biosynlhesis 
Ether lipid melabolism 
Sphingdipid malabolisfFi 
GlLJCO^inclate biosynlhei&is 
Indole alkaloid bip£ynihe$i$ 

One ca^an pool by folate 
Zeatin biosynthesis 
Pantothenato and CoA biosynlhesis 
Ubiqumono and other terporri>ld-quinone biosynthesis 
Valine, leucine and isoleucme biosynthesis 
Chloroalkane and chloroalkene degrBdalion 
Carbon lix alien pathways in prokaryotes 
Arachidoni£ acid metabolisin 
Stenrjid hormono biosynthesis 
TTiiamine metabolism 
Retinol melaboli$m 
^olenoocmpound melabolism 
Phgsphahdylinosilel signaling system 
Nipotinate and nicotinamide metabolism 
Glycosphiingol'jpiid biosynthesis -globo senes 
P^nicilliin and cephalosponn biosynthesis 
DOutpmine and D-glutamato metabolism 
Qelalain biosynthesis 
N-Qlycan biOsynth-eSlS 
Ciliate cycle (tCA cycle) 
Other glycan degradatiofi 
Polyketide sugar unit biosynthesis 
Photosynthesis 
Naphthalene degradation 

Tauinne and hypotaunne metdboii$in 
Foilafte biosynthesis 
Histiidin« metat>olisin 
Biosynthesis of snsamycms 
Upopoly saccharide biosynlhesis 
BiosynViesis of vancomycin group antibiotics 
Styiene degradation 
Tetracycline biosynthesis 
Steroid biosynthesis 
Bisphenol degradation 
Lysine biosynthesis 
Glycoiylplioisphalidy linos itoi«?P ijhsnchor bios y nthesis 
Cb-Branched dibasic acid melabo1i$ni 




10 



15 



20 



25 



30 



35 



40 



Number af sequences in pathways 

Figure 1 Pathways from the top ICQ up-regulated genes in TO interaction. Pathways were identified considering TO interaction. Amount of 
genes belonging to patliways is specified in eacli line. 
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KEGG pathways (down -regulated genes) 



SbRh ind »u£ros« m«laboiism 



Glyoxylala and die artioxy late melabaliain 
Thiamine mHlBboliann 
PlMfiylifllinln*, tynHliM and tryplochan bioftyntrmii 
InOBitel phosphatn nHlaboliflm 



pyruwate melabQtisni 



Gaia£l»*« iifk«iaboli »m 
QlycolyBis I Oucoinn>o«neBliB 
Nitn^en ndtivboiitni 
Aigrinlnt and pvi\%m mH^tKAtim 
PlMnylpfopanoid biosy nlheatn 
AmlncHMnzoaEa degradallorv 
DlXtg (ncKJttOJijfn . eytOCtwomt Pi^ 



AtccrbvM 4r>d aid*r*(« metaboiitni 

Fatty acid nnfltBt>olism 
jiyplcptian rnalabdicn 
Amln^ sugar ani3 nucl«Ollid«' ^U^ar melabolixri 
Meiaballam af xenobiollca by ^ytoclMome 

$laiicHd tKifmans E>lQ^ynlheai>s 
Arnin&acvURW*! tie»v"ihettt 
CvtHTn (VxHtion paltvways In prnKaryoteK 
FlavQnotff i>lQ»yn|rieskt 
Phanylalanina nialvboiktm 
CHlBte cycle (TCA cycle} 
{^UtBthiona msldbfillvm 
P«frto*« photpri*i« pdhMiy 



CartKMi flixatlon in plucyiasyntrialic aiganlAma 
Oxidative piKisphQryialion 
ctuKMHiitaiM and cmciraaiiiitnt d«ijrad4ii«n 

Fatty acid bioayntheaia 
G<yCftiopiK*ph«lipid rrielibotitm 
Indola aiKaloid biosyntheftia 
IJ|}lquln«tfl and oIlHr lerpanctd-qui nana bioaynthesliB 
Catf0<n« meiaboiisni 
Aracmidanic acid malvbotitnt 
Linoleic acid mfllabollBm 
Stynne deflradatlon 
Pttotptiaf Oylllnp^lliQl jl^nalinq t yil«im 
Glycine, aerlna and Ihracnine mHlat^oliam 
Oy&lB'nB and melFiiontne [ii>fltat>olism 
Lysine de^radiaiicn 
tHil»^lanirk« ni»1»boiiSrti 
Llmtmana and pinene dsgradaEion 
Geranial degiadalmn 
NaphlhaltM dA^rddalioni 
StUfur nuriabxriiam 
Histidine mdMralisni 
BloayntheaiB of ansamiyclna 
Ffucltraa and mannose nMlatwIiflm 
Riboflavin meiatpiit'n 
AlarMM, aipartalfr and glutamale meiatKilitm 
Selenacom pound [nHtat>oli am 
Aninocyanin biosyntneaia 
Pep<lid(>glycan blosynthesm 
tftDquinollnH alKalaid biaaynthesm 
$ptilngDllpld msljAqlisni 
Giyc»piilngoNpld biotyrvtttfltii) ^ giobo tenet 
TtopaiM, pIpMdirM afkd pyfldina alhalcM biuynth*4lft 
Panlclilln and ctpMlPtponn biPtynlhttit 
TaiTMneild b^chbonv bioay nthesrs 
Dtterpenwd bioByntheais 
Oain^encHd biosyninesis 
H-^Myttn frioivmnent 
Otttartypea of O^lycan bioaynthesia 
Pftmaiy blie acid t>iQ»yn|h99ii$ 
S-bsptpntyciA bietynlhetit 
Butiro^ln and neomycin NaiiyTith^ftia 
Valine, leucine »nd itdieudne biosynthesis 
Tpiuene dtgradiailen 
Biosynthesis -af unrsaturated latty aclda 
dycwamlneglycan tj^q^synthesls ■ heparan »uilale 
Fetalfr blesynlhosis 
Llpopotysac cttahda biosynthesis 
Flivone and rawonoi i>Josynthesis 
Bulanp«tt malabPtiSTi 
T cAll mceptor aignaling pa/thway 
ChlHocyclatieKane and chlaiotienzene degradallDn 
Fluoreb«rkZO«ta dagradtallcn 
Lysine blosynthesia 




0 5 10 15 20 25 30 

umtier of sequ«n«5 in paihwsy^ 

Figure 2 Pathways from the top ICQ down-regulated genes in TO interaction. Pathways were identified considering TO interaction. Amount 
of genes belonging to patliways is specified in eacli line. 
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Table 3 List of up-regulated genes observed for the TOR_TOS Interaction with respective annotation and 



fold-change values 


Name 


Blast Sol network genomics 


e-value 


Annotation^ 


Fold change (TOR X TOS) 


QUI 


SGN-E628893 


0 


— NA— 


1000,27 


0U2 


SGN-El 352064 


5e^^ 


Caffeine synthase 


642,35 


0U3 


SGN-El 326397 


0 


Acidic endocliitinase se2 


447,45 


0U4 


SGN-El 31 0344 


G 


Metallotliionein-lil<e protein 


280,91 


0U5 


SGN-El 334735 


G 


Kunitz trypsin inliibitor 


238,17 


0U6 


SGN-E1316291 


0 


Acid pliospliatase 


236,53 


0U7 


SGN-E661231 




Polygalacturonase-1 non-catalytic subunit beta 


236,37 


0U8 


SGN-El 32761 5 


0 


Kunitz trypsin inliibitor 


168,12 


0U9 


SGN-El 337775 


G 


Organ-specific protein 


155,92 


0U10 


SGN-E659257 


G 


Polygalacturonase-1 non-catalytic subunit beta 


152,17 


0U11 


SGN-El 352070 


G 


Caffeine synthase 


148,75 


0U12 


SGN-E642649 


0 


Cytokinin oxidase 


143,25 


GUIS 


SGN-E833713 




Asrl protein 


119,38 


0U14 


SGN-E667484 


G 


Protein 


115,77 


0U15 


SGN-El 31 9644 


G 


pr-10 type pathogenesis-related protein 


98,14 


0U16 


NM** 




Protein 


97,12 


0U17 


NM 




Protein 


94,57 


GUIS 


SGN-El 321 440 


G 


Acid phosphatase 


93,94 


GU19 


NM 




Invertase pectin methylesterase inhibitor family protein 


82,40 


GU2G 


SGN-E1312621 


G 


— NA— 


76,83 


GU21 


SGN-El 309331 


0 


mpbq msbq methyltransferase 2 


75,15 


GU22 


SGN-E837532 




Class iii chitinase 


75,05 


GU23 


SGN-E1316252 


0 


Heat shock 


73,55 


GU24 


SGN-El 325880 


0 


Swib complex baf60b domain-containing protein 


72,42 


GU25 


SGN-E832873 


le"^ 


Kunitz trypsin inhibitor 


63,26 


GU26 


SGN-El 3221 00 


le-^^ 


Polyphenol oxidase 


60,83 


GU27 


SGN-E660241 


G 


Polyphenol oxidase 


57,50 


GU28 


SGN-E671322 


0 


60s acidic ribosomal protein pO 


56,95 


GU29 


SGN-El 128614 


8e-" 


Protein 


53,42 


GU3G 


SGN-E682004 


0 


Lipid transfer protein 


50,27 


GU31 


SGN-El 334549 


G 


Type ii proteinase inhibitor family protein 


49,57 


GU32 


SGN-El 33771 5 


G 


Protein 


49,57 


GU33 


SGN-E640935 


0 


Class iii chitinase 


48,21 


GU34 


SGN-E990795 




Microsomal glutathione s- 


47,84 


GU35 


SGN-E668445 


0 


Serine-type endopeptidase inhibitor 


47,81 


GU36 


NM 




Metallothionein-like protein 


47,79 


GU37 


NM 




Polyphenol oxidase 


46,30 


GU38 


SGN-E640935 


G 


Class iii chitinase 


45,63 


GU39 


SGN-E657601 


G 


Protein 


44,71 


GU4G 


SGN-E835025 




4-hydroxy-3-methylbut-2-enyl diphosphate reductase 


44,21 


GU41 


SGN-El 3201 97 


G 


Peroxisomal membrane 


43,89 


GU42 


SGN-E636199 


G 


— NA— 


43,72 


GU43 


SGN-El 333755 


G 


Type ii proteinase inhibitor family protein 


41,88 
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Table 3 List of up-regulated genes observed for the TOR_TOS Interaction with respective annotation and 
fold-change values (Continued) 



0U44 


SGN-El 337775 


0 


Organ-specific protein 


41,68 


0U45 


SGN-E912118 


le-2^ 


Tartrate-resistant acid pliospliatase type 5 


41,04 


0U46 


SGN-E628829 


8e-^^ 


Cell wall protein 


40,77 


0U47 


SGN-E449176 




Phospholipid glycerol acyltransferase family protein 7 


39,90 


0U48 


SGN-E830846 


0 


r3h domain containing 


39,64 


0U49 


SGN-E628829 


8e^^ 


Oxygen-evolving enhancer protein chloroplast 


39,55 


0U50 


NM 




Auxin-independent growth promoter protein 


39,13 


0U51 


SGN-El 34931 2 


0 


Lipid transfer protein 


39,09 


0U52 


SGN-E639273 


^^-160 


Protein kinase domain containing expressed 


38,61 


0U53 


SGN-El 352095 


0 


Protein 


35,83 


0U54 


SGN-E838896 


0 


Cytochrome p450 


34,98 


0U55 


NM 




nadh dehydrogenase subunit f 


34,83 


0U56 


SGN-E838821 


0 


mta sah 


34,62 


0U57 


SGN-E669832 


0 


Cytokinin oxidase 


33,83 


0U58 


SGN-E643214 


0 


vesicle-associated membrane protein 714 


33,06 


0U59 


SGN-El 346029 


0 


Protein 


32,83 


0U60 


SGN-El 348577 


0 


Protein 


32,50 


0U61 


SGN-E674849 


le^^ 


Mitochondrial chaperonin hsp60 


32,40 


0U62 


SGN-E788243 


6e-2° 


Tartrate-resistant acid phosphatase type 5 


32,40 


0U63 


SGN-El 09388 


3e-^^ 


Protein 


31,93 


0U64 


SGN-El 326628 


0 


Formate dehydrogenase 


31,83 


0U65 


SGN-El 333903 


0 


Elongation factor-1 alpha 


31,16 


0U66 


NM 




Protein 


30,56 


0U67 


SGN-El 321 81 2 


0 


Peptidylprolyl isomerase 


30,27 


0U68 


SGN-El 31 0338 


0 


60s ribosomal protein 


30,00 


0U69 


NM 




Conserved hypothetical protein [Ricinus communis] 


29,87 


0U70 


SGN-E681444 


0 


Homeobox-leucine zipper protein 


29,71 


0U71 


SGN-El 168096 




Protein 


29,31 


0U72 


SGN-E818914 




Gibberellin 20 


29,19 


0U73 


SGN-E431715 


7e" 


Cellulose synthase 


28,99 


0U74 


SGN-E650445 


0 


Poly -binding protein 


28,29 


0U75 


SGN-E631106 


3e-^^ 


orf i polyprotein 


28,09 


0U76 


SGN-El 314273 


0 


Ankyrin repeat domain 


27,84 


0U77 


SGN-El 31 5499 




Atapy2 atpase nucleotide diphosphatase 


27,54 


0U78 


SGN-E1315958 


0 


nadh ubiquinone oxidoreductase bl4 subunit 


26,38 


0U79 


SGN-E830665 


0 


sec61 transport protein 


25,52 


0U80 


SGN-E1316141 


0 


Flavanone 3-hydroxylase-like protein 


24,98 


0U81 


NM 




cbl-interacting serine threonine-protein 


24,92 


0U82 


SGN-El 342733 


0 


Transcription factor lim 


24,70 


0U83 


SGN-E837532 


8e^° 


Class iii chitinase 


23,94 


0U84 


SGN-El 322208 


5e-9o 


Cysteine proteinase 


23,93 


0U85 


SGN-El 322866 




Protein 


23,85 


0U86 


SGN-El 31 5443 


0 


Acid phosphatase 


23,69 


0U87 


SGN-E662706 


0 


gdsl-motif lipase hydrolase family protein 


23,27 
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Table 3 List of up-regulated genes observed for the TOR_TOS Interaction with respective annotation and 
fold-change values (Continued) 



0U88 


SGN-E266690 


7e-" 


Protein 


23,16 


0U89 


SGN-E1319812 


0 


Protein 


22,87 


0U90 


SGN-El 334410 


0 


Protein 


22,82 


0U91 


SGN-E1321222 


0 


Triosepliospliate isomerase 


22,62 


0U92 


SGN-E648331 




Serine-tlireonine protein plant- 


22,61 


0U93 


SGN-El 351 186 


0 


Class iii chitinase 


22,36 


0U94 


SGN-E951741 


26^2 


Late embryogenesis abundant protein Ieal4- 


22,06 


0U95 


SGN-E686943 


0 


Dehydrin 


21,99 


0U96 


SGN-E674268 


le^° 


mta sah 


21,79 


0U97 


SGN-El 352075 


0 


7-methylxanthine n-methyltransferase 


21,73 


0U98 


SGN-El 323598 


0 


Transcription initiation factor lib 


21,29 


0U99 


SGN-El 348381 


le-^^^ 


Translation factor 


21,07 


OUlOO 


SGN-El 322408 


0 


Beta-glucosidase-like protein 


20,89 



*The annotation of each sequence was established on the Coffee Genome Database [31]. 
**No match with any sequence on the Solanaceae Genomic Database. 



fromlignin and isoflavones synthesis, downstream metabo- 
lites, such as cynnamyl alcoholdehydrogenase (-2,79) and 
isoflavone reductase (-1,20), are down-regulated in resist- 
ant plants. On the other hand, genes from biosynthesis of 
anthocyanins and tannins, such as flavonoide 3'-hydroxy- 
lase (TO = 24; Tl = 5; T2 = 8) and leucoanthocyanidin 
dioxygenase (TO = 2; Tl = 9; T2 = 3) are up-regulated at all 
times in resistant plants. This expression profile suggests 
that synthesis of anthocyanins and tannins is favored over 
synthesis of flavones. Also, activation of upstream genes 
such as phenylalanine ammonia lyase, chalcone synthase 
and flavonone dehydrogenase is observed at final stages 
in susceptible plants, indicating that phenylpropanoid 
biosynthesis is delayed. 

Several other defense-related genes are also positively 
regulated in resistant plants, including herbivore-response 
related genes glutathione-S-transferase and cysteine pro- 
teinase inhibitor, Apoptosis-related genes have a variable 
expression profile: catalase is up-regulated throughout 
insect development but citocrome c oxidase, superoxide 
dismutase and a senescence- associated protein have no 
differential expression, and polygalacturonase is up-reg- 
ulated only at TO. 



Validation of expression profile for selected 
candidate-genes 

We selected 18 genes for validation, listed with corre- 
sponding expression levels on Table 5. Genes exhibited a 
consistent expression pattern when quantified by either 
microarray or qPCR, and the Pearson coefficient for this 
comparison is 0.92 

Further qPCR analyses were performed to validate 
expression of selected leaf-miner resistance-associated 



candidate-genes (Figure 8). These included genes from 
pathways described above and genes with either no 
significant hits or similarity to unknown proteins, which 
may represent coffee specific genes, not yet identified or 
characterized. 

The putative caffeine synthase gene exhibited the greatest 
response to leaf-miner infection (Figure 8), as its expression 
was significantly higher (Relative Quantification value - 
RQv 230.45) in control resistant leaves, and also in later 
stages (RQv 155.24). The expression of gene SGN-E628893, 
encoding an unknown protein, is significantly higher (RQv 
1000) in control resistant leaves than in susceptible ones 
(Figure 8). However, transcript levels dropped in resistant 
plants at first stages of infection, such as oviposition 
(RQv -1.53) and egg ecclosion (RQv 1.22). This gene is 
a good candidate for differentiation of resistant and 
susceptible plants, although possibly is not related to 
resistance. 

The isocitrate lyase gene is down-regulated (RQv -6.38) 
in resistant leaves at TO. At initial steps of infection 
its transcript levels increased rapidly (RQv 5.49), but 
decreased during feeding stages. Other evaluated genes 
also exhibited a regulated expression. The glycerol-3- 
phosphate acyltransferase 6 gene, a trans-membrane 
protein associated with synthesis of cutin, is up regulated 
(RQv 11.12). A gene encoding metallothionein-like pro- 
tein, a class of metal-chelators proteins with possible 
anti-oxidant role, is also up-regulated (RQv 12.88) at 
initial stages of insect development. 

Genes encoding unknown proteins with lipase pro- 
tein domain, such as SGN-U585128, SL2.40ch08 and 
SGN-U585128 were activated in resistant plants, being 
up-regulated at oviposition and/or egg ecclosion with 
RQ values of 19.72, 19.03 and 10.30 respectively. As 
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Table 4 List of down-regulated genes observed for the TOR_TOS Interaction with respective annotation and 
fold-change values 



Narn6 


BIdst Sol n6twork Q6nornics 


6-valu6 


Annotation* 


Fold chdnQ6 
(TOR X TOS) 


0D1 


SGN-El 320843 


0 


Protein 


-445,87 


0D2 


SGN-E676870 


0 


Protein 


-235,45 


0D3 


SGN-El 320843 


0 


Zinc finger 


-126,65 


0D4 


SGN-El 325444 


0 


Hypotlietical protein VITISV_000181 [Vitis vinifera] 


-62,54 


0D5 


SGN-E835732 


0 


Tapetum-specific protein lla-115 


-49,63 


0D6 


SGN-E661762 


0 


— NA— 


-35,43 


0D7 


NM** 




PREDICTED: liypotlietical protein [Vitis vinifera] 


-24,73 


0D8 


NM 




Glycerol-3-pliospliate acyltransferase 6 


-17,93 


0D9 


SGN-El 321 887 


0 


gdsl-motif lipase hydrolase-like 


-16,98 


0D10 


SGN-El 033676 


5e-^^ 


er glycerol-phosphate acyltransferase 


-16,50 


0D11 


SGN-El 322050 


0 


Extensin-like protein 


-14,88 


0D12 


SGN-E791894 


le"^ 


Cytochrome p450 


-12,96 


0D13 


SGN-E837009 


0 


Cytochrome p450 


-11,39 


0D14 


NM 




Cytochrome b 


-11,12 


GDI 5 


SGN-El 349845 




gdsl-motif lipase hydrolase family protein 


-10,27 


0D16 


SGN-El 31 8049 




Gibberellin-regulated protein 1 


-9,89 


GDI 7 


SGN-E660879 


0 


Zinc finger 


-9,88 


GDIS 


SGN-E835266 


0 


Isocitrate lyase 


-9,22 


GDI 9 


SGN-E673783 


0 


Isocitrate lyase 


-8,84 


GD2G 


NM 




— NA— 


-8,83 


GD21 


SGN-E660879 


0 


Zinc finger 


-8,73 


GD22 


SGN-E680272 




Protein 


-8,498 


GD23 


SGN-E898278 




Protein 


-8,43 


GD24 


SGN-El 328871 


0 


— NA— 


-8,16 


GD25 


NM 




Cytochrome p450 


-8,02 


GD26 


SGN-E678498 


0 


Protein 


-7,36 


GD27 


SGN-El 31 9644 


0 


Serine-threonine protein plant- 


-7,35 


GD28 


SGN-E838812 




— NA— 


-7,34 


GD29 


NM 




Pathogenesis-related protein 1 


-7,27 


GD3G 


SGN-El 31 2882 


0 


Heat shock protein 


-7,17 


GD31 


NM 




— NA— 


-7,12 


GD32 


SGN-FR^DRDfi 


0 


r~\/tnrhrnrnp n4-Sn 


-7 06 




J^lN II 1 JJ^UUZ 


n 
u 


dUL LI dl IbpUl Lcl 


-7 nn 
/ ,uu 


GD34 


SGN-F6S9349 


0 


f-ili itrithinnp '^-trsn'^fprs'^p n<^ti \f) 

vJIULULI ll^l IC: D Liai IJlCTiajC: yjLUVJ 


-6,85 


GD35 


buN-L 1 3 1 / 1 04 


0 


Aspartyl protease family protein 


-6,74 


GD36 


SGN-El 322588 




atlg72120 f28p5_2 


-6,53 


GD37 


NM 




— NA— 


-6,46 


GD38 


SGN-El 350292 


0 


Lactoylglutathione lyase family protein 


-6,46 


GD39 


NM 




Achain crystal structure of a cell-wall invertase from 
Arabidopsis thaliana in complex with sucrose 


-6,36 


GD4G 


SGN-El 335955 




Retroelement pol polyprotein 


-6,24 


GD41 


SGN-E820310 




Xyloglucan endotransglucosylase hydrolase protein 22 


-6,23 


GD42 


SGN-E836814 




Leucine-rich plant specific 


-6,20 
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Table 4 List of down-regulated genes observed for the TOR_TOS Interaction with respective annotation and 
fold-change values (Continued) 



0D43 


SGN-E686810 


0 


Zinc finger 


-6,09 


0D44 


SGN-E839045 


0 


Glucose-metlianol-clioline oxidoreductase family protein 


-6,06 


0D45 


SGN-E531670 


le-^2 


Hydrolyzing o-glycosyl 


-5,97 


0D46 


SGN-E1216540 




Aminotransferase family protein 


-5,91 


0D47 


SGN-E747084 




Alkaline alpha-galactosidase seed imbibition protein 


-5,81 


0D48 


SGN-E836814 




Leucine-rich plant specific 


-5,71 


0D49 


SGN-El 345225 


0 


— NA— 


-5,71 


0D50 


SGN-El 325272 


0 


Protein 


-5,69 


0D51 


SGN-E626178 






-5,62 


0D52 


NM 




Alkaline alpha-galactosidase seed imbibition protein 


-5,60 


0D53 


NM 




CinnamoyI reductase-like protein 


-5,57 


0D54 


SGN-El 349228 


0 


Proline dehydrogenase 


-5,54 


0D55 


SGN-E775239 




Kinesin like protein 


-5,46 


0D56 


NM 




Methyl-accepting chemotaxis sensory transducer 


-5,45 


0D57 


NM 




Tr?in<^rrintinn fartnr 


-5,42 


0D58 


NM 




D 1 lyu 1 UAy loUUU Ly 1 d Lc Ucl lyuiuytrl Idoc Idlllliy piULclll 


-5 39 




NM 




wULcl 1 1 Icl 1 lUI dl Ic puilll piULclll 




0D60 


J'OIN l_ 1 J 1 / O-JJ 


Q 


1 1 dl loLI IjJLIUI 1 dLLIVdLUI 


-5 31 


0D61 


SGN-E666413 


0 


bahd family clade v 


-5,20 


0D62 




0 


r^|\/r-prnl-^-nhnc;nh;^tp Hph\/H rnnpn;^<^p 


-5 18 


0D63 




0 


1 Ul lUjJIdoL IIILIIIIolL 


-5 1 7 


0D64 


SGN-1 ^^1467 


0 


Arh;^in rf\/'^t^\ c;triirtiirp nf ^ rpll-\A/;^ll in\/prt;^c;p frnm 

arabidopsis thaliana in complex with sucrose 


-5 08 


0D65 


NM 




Inner-membrane translocator 


-5,00 


0D66 


NM 




Beta-ig-h3 fasciclin 


-4,94 


0D67 


SGN-E653486 




bahd family clade v 


-4,90 


0D68 


NM 




Stachyose synthase 


-4,86 


0D69 


SGN-El 3491 01 




60s acidic ribosomal protein pi 


-4,85 


0D70 


NM 




Disease resistance 


-4,79 


0D71 


SGN-E666413 


0 


bahd family clade v 


-4,75 


0D72 


SGN-El 31 3854 


0 


Cytosolic aldehyde dehydrogenase 


-4,75 


0D73 


SGN-El 320568 


0 


Protein 


-4,72 


0D74 


SGN-E748200 


le-22 


Anthranilate synthase alpha subunit 


-4,71 


0D75 


SGN-El 31 6428 


0 


Heat shock protein 


-4,68 


0D76 


NM 




Protein 


-4,67 


0D77 


SGN-E700960 


8e-^^ 


Magnesium transporter 


-4,67 


0D78 


SGN-El 321 133 


0 


Glutathione s-transferase 


-4,66 


0D79 


SGN-El 309644 


le^^ 


— NA— 


-4,66 


0D80 


NM 




nac domain ipr003441 


-4,63 


0D81 


SGN-E791702 


0 


Zinc finger 


-4,57 


0D82 


SGN-E524668 


0 


Protein 


-4,49 


0D83 


SGN-E667829 


0 


gdsl-motif lipase hydrolase family protein 


-4,46 


0D84 


SGN-E955597 


4e^^ 


Proline dehydrogenase 


-4,44 


0D85 


NM 




— NA— 


-4,43 
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Table 4 List of down-regulated genes observed for the TOR_TOS Interaction with respective annotation and 
fold-change values (Continued) 



0D86 


SGN-E1312314 


0 


PREDICTED: hypothetical protein [Vitis vinifera] 


-4,39 


0D87 


NM 




Undecaprenyl pyrophosphate phosphatase 


-4,37 


0D88 


SGN-El 350610 


0 


Cytosolic class i small heat-shock protein 


-4,35 


0D89 


NM 




— NA— 


-4,34 


0D90 


SGN-E747084 




Alkaline alpha-galactosidase seed imbibition protein 


-4,34 


0D91 


NM 




Protein 


-4,32 


0D92 


SGN-E528554 




pili assembly chaperone 


-4,31 


0D93 


SGN-E528554 


7e-3o 


Cell-wall invertase 


-4,28 


0D94 


NM 




nbs-Irr resistance protein 


-4,26 


0D95 


NM 




Transcriptional family 


-4,26 


0D96 


NM 




— NA— 


-4,24 


0D97 


NM 




Oxysterol binding protein 


-4,21 


0D98 


SGN-El 196563 


2e^^ 


ap2 domain-containing transcription factor 


-4,20 


0D99 


SGN-E834183 


0 


nac domain protein nac2 


-4,19 


ODIOO 


SGN-E678677 


0 


— NA— 


-4,19 



*The annotation of each sequence was established on the Coffee Genome Database (Vieira et a!., 2006). 
*^No match with any sequence on the Solanacea Genomic Database. 



these genes are down-regulated on leaves at TO, this 
profile during insect development suggests a role during 
defense response (Figure 8). 

The expression profile of polygalacturonase gene was 
not confirmed by real-time PGR. According to microarray 
in silico analysis this gene is activated at TO (fold change 
value 236,370), and afterwards transcript levels drop 



along insect development (fold change value 1,628). Yet 
by real-time PGR analysis transcript levels are similar 
in susceptible and resistant plants at all times. This result 
may reflect differences on genetic background of evaluated 
plants, once they are part of a population still segregating 
for some characteristics. However, as this was the only 
observed discrepancy between all performed analyses. 
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Figure 3 Biological process for up-regulated genes. Functional categorization of up regulated genes (%) in coffee genotypes (susceptible and 
resistant to leaf miner) without insect. 
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Figure 4 Biological process for down regulated genes. Functional categorization of down-regulated genes (%) in coffee genotypes (susceptible 
and resistant to leaf miner) without insect. 



these genetic differences may not be associated with the 
resistance trait. 

Discussion 

The use of resistant or tolerant cultivars represents an 
alternative for reducing the use of chemical defensives, 
the costs of production, and the negative impact over 
environment. In coffee, breeding for herbivore resistance 
is important once Coffea arabica, the main commercial 
Cojfea species, is susceptible to almost all known coffee 
pathogens. Then resistance must be transferred from 
other compatible species, which is not always possible 
due to limited efficiency of inter-specific crosses. 

At the Agronomic Institute (Campinas, Brazil) leaf-miner 
resistance genes have been transferred from C. racemosa 
to C. arabica through traditional breeding strategies, and 
although the program is currently at advanced generations, 
no resistant commercial cultivar is available yet. The lack of 
knowledge of molecular aspects controlling the resistance 
response, and the restricted genetic variability of breeding 
populations [34] limit the development of genomic-based 
selection tools. In this context, we aimed to provide 
information on molecular aspects of leaf-miner defense 
mechanisms and identification of reliable candidate-marker 
genes for assisted-selection. Those genomic tools associated 
with traditional breeding strategies guarantee that agrono- 
mical traits such as productivity and cup quality will 



be selected in advanced generations. Also, once novel 
genomic-markers associated with other desirable traits are 
developed for coffee genotypes, a genome-wide selection 
strategy will be possible to develop multiple-pathogen 
resistant cultivars. We chose the microarray analysis 
approach to assess the co-expression of a large amount 
of genes, including those that are not looked at in common 
analyses. Nevertheless, the results described in this work 
indicated that pathways regularly activated in response to 
herbivory, such as linoleic acid cycle, phenylpropanoids 
synthesis and apoptosis, are also activated during coffee- 
response to leaf-miner. Genes associated with jasmonate 
(JA) synthesis, such as lipoxygenase and enoyl-CoA hydra- 
tase, and with flavonoids synthesis, such as chalcone 
synthase and flavanone 3-hydroxylase-like, are up-regulated 
in resistant plants. Also, pathways from the primary me- 
tabolism, such as the citric acid cycle are down-regulated 
during leaf-miner defense response in resistant plants, a 
profile observed also in conifers [35]. 

Results of this study provide evidence that most genes 
encoding enzymes from the citric acid cycle are down- 
regulated in resistant plants (Figure 5). In a parallel ana- 
lysis, a metabolite profile was established for resistant and 
susceptible genotypes using an NMR-based technique 
[33] and indicated that malate levels on resistant leaves 
are lower than in susceptible ones. Malate results from 
conversion of either fumarate or glyoxylate. Expression 
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Figure 5 Citric acid pathway (KEGG database). Enzymes are identified by E.C. number and tliose corresponding to differentially expressed genes 
are highlighted according to the color scale. Color codes for each gene at TO are as follows: red for up-regulated, green for down-regulated and grey 
for no differential expression. E.C number and correspondent enzyme: 1.1.1.37 - Malate dehydrogenase; 1.1.1.41 - Isocitrate dehydrogenase (NAD 
(+); 1.1.1.42 -Isocitrate dehydrogenase (NADP(+); 1.2.4.1 - Pyruvate dehydrogenase (acetyl-transferring) (Pyruvate dehydrogenase) 1.3.99.1 - Succinate 
dehydrogenase (Fumarate reductase); 2.3.3.1 - Citrate (Si)-synthase, 4.1.1.32 - Phosphopyruvate carboxylase; 4.2.1.2 - Fumarate hydratase (Fumarase). 



of fumarase, that converts fumarate into malate, is down- 
regulated at TO in resistant genotypes (-2), what could 
explain the low malate levels. Production of malate 
from glyoxylate may also be deficient in resistant plants, 
once genes encoding for malate synthase, that converts 
glyoxylate into malate, and for isocitrate liase, the upstream 
enzyme that converts isocitrate into glyoxylate and succin- 
ate, are both down-regulated (-2,51 and -9, respectively). In 
contrast to this profile, susceptible plants exhibit a regular 
expression levels for these genes at TO. Therefore, both 
metabolic and transcriptional profiles support the affirma- 
tions that citrate cycle is down regulated in leaf-miner 
resistant coffee plants, and the model of down-regulation 
of primary metabolism in herbivore-resistant plants [36]. 

Biosynthesis of JA starts with alpha-linoleic acid release 
in non-injured tissues, triggered by systemin and phospho- 
lipase A2, Alpha-linoleic acid is then converted to JA after 
enzymatic steps performed by 13-lipoxygenase (LOX), 
allene oxide synthase (AOS), jasmonate o-methyltrans- 
f erase and others [37,38]. Several genes from the JA 



biosynthesis pathway are up-regulated in resistant 
plants at TO, including those from downstream steps 
such as jasmonate o-methyltransferase which expression is 
10-fold higher than in susceptible plants. All genes of the 
JA biosynthesis are either down-regulated or up-regulated 
at later stages in susceptible plants, as for instance expres- 
sion of LOX (TO -8.66; Tl 1; T2 2) increases only at Tl. 
These observations suggest that the JA signalling pathway, 
including intermediate signaling- molecules such as oxo- 
pentenyl-cyclopentane (OPC), may be impaired in suscep- 
tible plants. Down-regulation of genes from later steps of 
JA biosynthesis, such as allene oxide cyclase, allene oxide 
synthase, carboxyl methyltransferase, the enzyme that 
converts jasmonic acid into methyl jasmonate, is observed 
at Tl and T2 in resistant plants. Therefore, a feedback 
regulation may be activated, with a re-programming of 
transcriptional response upon leaf-miner infection. 

Genes associated with biosynthesis of secondary com- 
pounds are shown here to be regulated. Expression profile 
of genes from phenylpropanoids biosynthesis, both 



Cardoso et al. BMC Genomics 2014, 15:66 
http://www.bionnedcentral.conn/1 471 -21 64/1 5/66 



Page 14 of 20 



a-LINOLENIC ACID METABOLISM 



-I 



Range of Fold Change in the Pathway 



30-39 



Phosphat idy Ichol i ne 
311 32| 



a-Linolenic acid O ( FADS2 — Stearldonic acid 



Volicitin 



l7~Hydroxy~ 
Itnolenic acid 



2(R)-HpOTrEr 
Oi — <>^H DOX 1 



2(R)-H0TrE 



spontaneous 
Heptadecatrienal < 



ntan 

i 



6 

(8Z,ilZJ4Z)- 
Heptadecatrienoic acid 



9(S)-H0TrE 

2 



9(S)-HpQTrE 



0 3,6- 



3,6-Nonadjenal 



HPLl J-^^^—^O 9-Oxononanoic acid 
CoJneJenic acid 



9,10-EOTrEi 



-♦O lO-OPDA 



l3(S)-H0TrE 



HpOTrE 



3-Hexenal 

2 

12-Oxo-9<Z)- 
dodece noic acid 
HPU V —^O^ 5 3 3 - 



20*29 

2-9 

-2-9 
-10-19 

non differential 



12»13-EOTrE 



42.1.921 



*o 

Etherolenic acid 



^ 

'Traumatin Traumatic acid 



S3S9£ 



0PC8 



OPCLl 



OPCS-CoA 



J2-0PDA 



0PC4-C0A 



CX- H23.1.16} 
AGX I 



3-0x0- 
OPC6-C0A 



trans-2-Enoyl- 
OPC6-C0A 



I ACX 



(^-Oxidation 



-CX H23J.16I MFP2 I — Q 

OPC6-C0A ^-Oy^o- trans-2-EnDyl 



(+)^7-Iso 
Ja!»nnonate 



trans-2-Enoyl- 
OPC4-C0A 



j23.1.ig| — »Q 1 312.- I — Kp* 

3-Oxo- JA-CoA 
OPC4-C0A 



OPC8-C0A 
(-)-Jasmonate 

|21JJ4l|. 



OPCS-CoA 



om — 

(+)-7-IsoiTiethyl- 
jasmonate 



(-)-Methyl- 
Jasmonate 



Figure 6 Alfa-Llnolenic acid metabolism pathway (KEGG database). Enzymes are identified by E.C. number and tliose corresponding to 
differentially expressed genes are highlighted according to the color scale. Color codes for each gene at TO are as follows: red for up-regulated, green 
for down-regulated and grey for no differential expression. E.C number and correspondent enzyme: 1.3.1.42 - 12-oxophytodienoate reductase; 
3.1.1.4 - Phospholipase A2; 1.13.11.12 - Linoleate 13S-lipoxygenase; 5.3.99.6 - Allene oxide cyclase; 2.1.1.141 - Jasmonate 0-methyltransferase; 
FADS2 - fatty acid desaturase 2 ; IVIFP2 - 3-hydroxyacyl-CoA dehydrogenase/ enoyl-CoA hydratase. 



up-stream genes such as PAL and CHS, and downstream 
genes such as flavonoide S-hydroxylase, leucoanthocyanidin 
dioxygenase, reveals a preferential synthesis of tannins and 
anthocyanins instead of ligninin, flavones and isoflavones 
in resistant plants. This profile indicates a direct defense 
strategy against leaf-miner, once among anthocyanins and 



tannins are found toxic compounds with antifeedants 
effects over insects [39,40]. 

Another gene linked to secondary metabolism is a 
putative caffeine synthase, which encodes one enzyme from 
caffeine biosynthetic pathway. Caffeine is an alkaloid 
distributed in coffee plant tissues and organs. The fact 
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Figure 7 Phenylpropanoid biosysthesis pathway (KEGG database). Enzymes are identified by E.C. number and tliose corresponding to 
differentially expressed genes are highlighted according to the color scale. Color codes for each gene at TO are as follows: red for up-regulated, 
green for down-regulated and grey for no differential expression. E.C number and correspondent enzyme: 1.1.1.195 - Cinnamyl-alcohol 
dehydrogenase; 1.2.1.44 - Cinnamoyl-CoA reductase; 1.11.1.7 - Peroxidase; 1.14.13 monoxygenase; 2.1.1.68 Caffeate 0-methyltransferase; 
2.3.1.91 - Sinapoylglucose-choline 0-sinapoyltransferase 2.3.1.99 - Quinate hydroxycinnamoyltransferase, 2.3.1.133 - Shikimate 
0-hydroxycinnamoyltransferase; 4.3.1.24 - Phenylalanine ammonia lyase. 



that expression of a gene from its biosynthesis is significantly 
increased upon leaf-miner infection suggests that caffeine 
may have a role in defense response. However, several 
studies regarding caffeine and leaf-miner development 
indicated that this compound has no effect on insect 
survival rates [41,42]. 

Once the ultimate goal of this study is to identify poten- 
tial candidate for markers, several genes were selected 
for validation using real-time PGR. Potential candidates 
include: isocitrate lyase, which increased expression during 
initial steps of leaf-miner infection may be co-related with 
reduction of primary carbon metabolism; putative caffeine 
synthase, part of an important pathway of coffee plants; 
glycerol-3-phosphate acyltransferase 6, a gene associated 
with lipid metabolism and part of cutin biosynthetic path- 
way, a secondary metabolite [43]; and finally metallothio- 
nein-like protein gene, a metal-transporter protein family 
with an uncertain role in plant metabolism but previ- 
ously associated with redox responses [44]. Future analyses 
include cloning and re-sequencing genomic regions of 
target genes from different genotypes in order to identify 
suitable polymorphisms. 



Among selected genes are those that have no similarity 
with any known reported gene or protein. Although they 
could not yet be associated with a biological process, 
their expression profile was very specific and related to 
defense response. For instance, genes SGN-E676870, SGN- 
E1128614 SGN-E1320843 were activated upon leaf-miner 
infection in resistant plants, and therefore represent good 
candidates for further investigation. Another interesting un- 
known gene is SGN-E628893, which is highly activated in 
resistant plants at TO but is repressed upon infection. The 
expression profile indicates that this gene is useful for early 
differentiation between resistant and susceptible plants. 

In summary, differential expression profiles between 
resistant and susceptible genotypes are observed even in 
the absence of leaf-miner, indicating that defense is already 
build up in resistant plants, as a priming mechanism. 
Then, a systemic defense response may be more rapidly 
activated in resistant plants, once basic compounds such 
as nitrogen and sugars are readily available as a result 
of repression of primary metabolism. This shift in plant 
metabolism is common after a pathogen attack, where 
defense-related pathways are activated, resulting in 
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Table 5 Expression Pattern of eighteen genes selected for validation 



Name 


Blast SOL network 
genomics 


Annotation 


Fold-change value* 
(microarray) 


AACt value** Coef. Pearson.*** 
(relative expression) 


m 1 1 "7 
UU 1 / 


IN IVI 


Protein 


94 


u,yz4z 


0U6 


bGN t\5\ Dzy 1 


Acid pliosphatase 


236 


1 5 


ni 1 1 /I 
UU 14 


oUN Ldd/4o4 


Protein 


lie 


3D 


uuz 


j>oIn L I jjZUD^ 


Itrll It: byi ILl Idbc 


O^Z 


1 /y 


0U4 


buN b 1 0 1 (Jo44 


Metallothionein-lil<e protein 


280 


8,45 


UUozz 


oUN L4jUzz I 


nadp-dependent d-sorbitol-6-pliosphate 
deliydrogenase 


D,J0 


5 


QUI 


SGN E628893 


No hits 


1000 


186,9 


0U42 


SGN E636199 


— NA— 


43 


8,16 


0D6 


SGN E661762 


— NA— 


-35 


-4,3 


0D10 


SGN El 033676 


Glycerol-pliospliate acyltransferase 


-16 


-1,8 


0D2 


SGN E676870 


Protein 


-235 


-120 


GDI 


SGN El 320843 


Protein 


-445 


-200 


GDI 9 


SGN E673783 


Isocitrate lyase 


-8 


-6 


GD3 


SGN El 320843 


Zinc finger 


-126 


-17 


GD4 


SGN El 325444 


Hypotlietical protein 


-62 


-36 


GD8 


NM 


Glycerol-3-pliospliate acyltransferase 6 


-17 


-2,3 


GD9 


SGN El 321 887 


gdsl-motif lipase hydrolase-like 


-16 


-2,2 


GD7 


NM 


Hypothetical protein 


-24 


-5,8 



*Microarray analysis value. 
*^qRT-PCR analysis value. 
*^*Pearson coefficient correlation value. 



reduction of growth and reproduction, and in changes 
on Unk-source relationship [45]. During herbivore- 
defense transcript levels of genes involved in photosyn- 
thesis are also down-regulated [21,46], probably as a 
strategy to liberate nitrogen compounds for the secondary 
metabolism. Maintenance of these physiological and meta- 
bolic states has a high energetic cost, and could represent 
a survival limitation if nutritional conditions on the 
field are depleted. Actually, field observations in culti- 
vated areas demonstrated that leaf-miner resistant coffee 
plants, under a severe nutritional deficit, are attacked by 
the leaf-miner at the same intensity as susceptible plants 
[27]. 

Conclusions 

As a basal defense state is decisive for triggering a rapid 
resistance response, genes associated with priming validated 
here, represent key genes for assisted-selection. Future 
studies will focus on comparisons of selected genes gen- 
omic sequences, from both resistant and susceptible paren- 
tal lines, to identify suitable marker polymorphisms. 

Methods 

Plant materials 

Resistant and susceptible coffee progenies were developed 
by the Coffee Breeding Program from the Agronomic 



Institute (I AC), Campinas, Sao Paulo, Brazil. The evalu- 
ated population (H14954-46), with 136 plants, is a F2BC5 
generation of the inter-specific cross (C. racemosa X C. 
arabica) X C. arabica. 

Plants were evaluated regarding the defense- response 
to Leucoptera coffeella using infestation methodology 
described by Guerreiro-Filho et al. [28]. After egg hatch- 
ing, 1.8 cm leaf discs were taken from the leaves using 
a cork bore. Discs were placed on damp plastic foam 
and maintained in a plastic box for two weeks. Resistance/ 
susceptibility response was visually scored according to 
the evaluation scale defined by Ramiro et al. [26]. 

Coffee leaf miner infestation of selected plants 

Fifteen resistant (R) and fifteen susceptible (S) coffee plants 
previously selected were used for leaf-miner infestation. 
Plants/seedlings of each group (R and S) were challenged 
with L coffeella in rearing cages and following the same 
procedures described above. Control non-infestaded plants 
of each group were also evaluated. Three independent 
infestation experiments were used for further analyses. 

Leaves were collected from the third and fourth pair 
from plants during different stages of insect development. 
Stages corresponded to egg hatching, after 1 to 5 days 
after infestation with L, coffeella (Tl), and egg eclosion 
and tissue injury, after 6 to 10 days after infestation (T2). 
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Correiation between foldchange values 
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change values from microarray analysis and soft color from qRT-PCR. 



Control non-infested leaves (TO) were also collected for 
each genotype. Three leaves of each plant and each stage 
of insect development were collected (totalizing nine leaf 
per time of sampling/genotype) and immediately frozen in 
liquid nitrogen and stored at -80°C until RNA extraction. 
Experimental design was completely randomized includ- 
ing three replicates for each sample. 



RNA isolation and preparations 

Total RNA for both NimbleGen microarray hybridization 
and real-time qPCR experiments was isolated using 
protocol described by Chang et al. [47]. RNA extractions 
were performed using 2 g of tissue of pooled samples. 
All RNA samples were analyzed by formaldehyde-agarose 
gel electrophoresis and by spectrophotometry to assess 
physical and chemical integrity. To avoid contamination 
by polyphenols, carbohydrates and proteins, only RNA 
samples with OD 260/280 and 260/230 > 1.8 were selected 
for further analysis. For microarray hybridizations, extracted 
RNA was also checked for purity and degradation using an 
Agilent Bioanalyzer 1000 (Agilent Technologies). Samples 
were stored at -80°C until further use. 



cDNA double strand synthesis, labeling and hybridization 

Ten thousands nanograms (10.000 ng) of each RNA 
sample were pooled and treated with DNAse- RNAase 
free for cDNA synthesis and labelling. Three biological 
replicates of each treatment were used for hybridization 
with the cDNA microarray chip. Equal amounts of each 
replicate from resistant and susceptible plants were pooled 
respectively to minimize variation between individual 
RNA samples. All RNA samples were sent to Roche 
NimbleGen Systems, where cDNA synthesis and Cy3 
labeling were performed following the manufacturers 
procedures (Nimblegen Gene Expression Analysis proto- 
col, Nimblegen Systems, Inc., Madison, WI, USA). Equal 
quantities of total RNA of each sample were converted 
to double strand cDNA (cDNA Synthesis System, Roche 
Applied Science). All the required equipments, reagents 
and procedures were provided and executed by Roche/ 
NimbleGen. 

Design and production of the Coffea ssp. Nimblegen® 
custom array 

Arrays were designed using sequence information available 
at the Brazilian Coffee Genome Project, which contains 
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sequences of around 33 K genes identified in EST libraries 
prepared from different physiological and metabolic 
situations [31]. The Coffea dataset was composed by 
quality-filtered contigs from different non-normalized 
ESTs/cDNA libraries of two coffee species Coffea 
arabica, Coffea canephora and Coffea racemosa, and 
by singlets of this assembly. Only sequences with at 
least one blast hit against NR database (e-value <le-10) 
were used as source sequences to generate probes for 
the 12 coffee microarray. The probes were designed by 
Roche-NimbleGen software, which selected unique 
sequences regions for each gene to avoid multiple 
hybridization with gene family members. Each micro- 
arrays consisted of 135.000 probes with length of 48 
nucleotides and Tm average from 68°C to 76°C, represent- 
ing 22,000 genes, with a minimum of 6 probes/gene. The 
final probe list was submitted to Roche-NimbleGen, Inc. 
(Madison, WI, USA) for quality control and subsequent 
probe array layout. Additional probes were also included 
on the microarray by Roche-NimbleGen, Inc. for quality 
control of the hybridization process. Microarray 



manufacture was synthesized in situ by photolithography 
on glass slides using a random positional pattern by Nim- 
bleGen (http://www.nimblegen.com/). 



Normalization and statistical analysis 

Hybridized-microarray slides were imaged with a high 
resolution array scanner (GenePix 4000B Microarray 
Scanner, Molecular Devices Corp., Sunnyvale, CA, USA) 
and fluorescent signal intensities from each spot were 
quantified using NimbleScan Software (NimbleGen Sys- 
tems Inc.). The intensity values were normalized using 
the oligo package from R statistics software [48]. The 
workflow used to normalize our data was followed as 
explained by the package provided for Nimblegen® expres- 
sion microarrays. Fold change values were calculated 
comparing resistant and susceptible genotypes with and 
without infestation. All clusters were annotated using 
the blast2go software [48] in order to label them with 
their probable molecular function, biological process 
and cellular component. An automatic pipeline using Perl 



Table 6 Primer sequences used for validation in qRT-PCR analysis 



Name 


Blast SOL 
network genomics 


Annotation 


Foward sequence 


Reverse sequence 


0U17 


NM 




Protein 


ACTACCAACATOACAGCAGCTC 


™accctgtoaagg™gtgc 


0U6 


SGN 


E1316291 


Acid pliospliatase 


CTAATOACCCTCTCCGCATGAT 


GCCAACTCAGGCAATOTATACG 


0U14 


SGN 


E667484 


Protein 


TAGTCAAGAATATGGGCATGGAC 


ATACCTOTOATOACGCCTO 


0U2 


SGN 


El 352064 


Caffeine syntliase 


AAAGGGAGCAmACTCTOCAAAG 


AGCATGCATCCTGAGAAATGTGGTA 


0U4 


SGN 


El 31 0344 


Metallotliionein-lil<e protein 


ATOGTCTGCTCTGTGAAGATGT 


ATACATGmCCGCAGmCCT 


0U322 


SGN 


E450221 


nadp-dependent d-sorbitol-6-pliospliate 
deliydrogenase 


CCmGTGGCTOTAAGCAAAT 


GGAAAGCAGAGATOACAAACAG 


OUl 


SGN 


E628893 


No hits 


CAAGGAAGATGC^GACGAT 


TGTAATOTGCTGCTGGTGCTAC 


OUl 


SGN 


E628893 


No liits 


CAmAGmCGAAGGGGACAA 


GGATACAGCCGGTAGGACTAACT 


0U42 


SGN 


E636199 


— NA— 


ACCCGCCGGGAAACC 


GATGCACAGACAGGAATCACAAC 


0D6 


SGN 


E661762 


— NA— 


TOGTAATGATGGAAGTGTCCTC 


GCCAGTAATGGGATOTAGTOA 


ODIO 


SGN 


El 033676 


Glycerol-pliospliate acyltransferase 


CTAGCTOACCAGGAAAGACAAG 


GACTCAGGACTGCTCAmCA^ 


0D2 


SGN 


E676870 


Protein 


GGCAACTACTGCATOTATCAGC 


AAATGGATGAGCTGAAGGAGAAC 


ODl 


SGN 


El 320843 


Protein 


ACTAGTACTGGGTGTOCCTCAA 


GGTGAGCAAATAGTOTOTOC 


0D19 


SGN 


E673783 


Isocitrate lyase 


GGCCAGGAGCAACAGACA^ 


ATCTCTCACAATCTOACmCCA 


0D3 


SGN 


El 320843 


Zinc finger 


CTGATOCGTCCGCTATCTCA^ 


AACCTATCGGACCTGTACCTG^ 


0D4 


SGN 


El 325444 


Hypotlietical protein 


TACTGGCACTAATGGAGGAAATG 


AATGGGGACAGATGTATCATCAC 


0D8 


NM 




Glycerol-3-pliospliate acyltransferase 6 


TOAAGAGmGGTACTGACGTG 


ATCATGGTCTGTCTCTCGATCTC 


0D9 


SGN 


El 321 887 


gdsl-motif lipase hydrolase-like 


TCTOACTGGACTOCTCCAATG 


CCATOAAGmAGAGCCACAAC 


0D7 


NM 




Hypothetical protein 


CCAAGACAGTOATCTCCCTCTA 


CGTAGTAGCTAGATGGTGCCAGT 


0U9 


SGN- 


-El 337775 


Organ-specific protein 


GGmCmAGGGmCCTOCT 


CACAGIGIGIGIGI 1 1 IGI ICCT 


OUll 


SGN- 


-El 352070 


Caffeine synthase 


CCTAGCAAGCCA^GGAG 


ATOTOGCAAACCTGTGGA 


0U6 


SGN- 


-E1316291 


Acid phosphatase 


a™cggctatggcagaa™gc 


CACCATCmCCTOmGAGA 


0U3 


SGN- 


-El 326397 


Acidic endochitinase se2 


CAGCAAATOTOCCTATGTCC 


CAGCGmCAGGGTOACATAAG 


0U8 


SGN- 


-El 32761 5 


Kunitz trypsin inhibitor 


CTCTOCmCAmCTGCTCTO 


GACGTAGTACTCGACACCAGGAC 
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scripts was created to map each probe to its corresponding 
gene and annotation. 

Differentially expressed genes (fold change values between 
2 and -2) were identified using linear models and by 
taking into account technical and biological replicates. 
When individual probes met the criteria that average 
signals from resistant versus susceptible genotypes differed 
significantly by at least two fold, probes were selected for 
final analysis. 

Functional characterization of differentially expressed 
genes was performed using Blast2GO [49] and also through 
directed searches on Gene Ontology (www.geneontology. 
org), KEGG (www.genome.jp/kegg) e InterPro (http://www. 
ebi.ac.uk/interpro/) databases. 

Microarray validation 

Validation of selected differentially expressed genes was 
performed by real-time PGR. Gene-specific primers were 
designed using Primer Express 3.0 (Applied Biosystems) 
and Premier Primer 5.0 (Premier Biosoft International, 
Palo Alto, CA, USA). Gene sequences were aligned with 
GeneBank reference sequences using the tBLASTx tool. 
Possible ORFs and functional and conserved domains 
were identified using the Open Reading Frame Finder 
(ORF FINDER) and CDD tools from the NCBI database. 
In order to guarantee gene-specificity and avoid amplifi- 
cation of multigene families, primers were designed upon 
target regions which included the conserved domain and/ 
or motif and anchoring outside the conserved region. A 
list of designed primers is shown on Table 6. 

Corresponding cDNAs were synthesized from 250 ng 
of total RNA using the RevertAidTM Minus First Strand 
cDNA Synthesis Kit (Fermentas) according to the manu- 
facturer protocol. PGR products were amplified using 
primers designed with Primer Express 3.0 (Applied Bio- 
systems) and analyzed by Premier Primer 5.0 software 
(Premier Biosoft International, Palo Alto, CA, USA). In 
order to confirm primer specificity and presence of single 
amplicons, all PGR products were analyzed through a 
dissociation curve, with temperature varying from 60°C 
to 95°C. 

Thermocycling and fluorescence detection were per- 
formed using ABI Prism 7300 Sequence Detection System 
(Perkin-Elmer Applied Biosystem). Real-time PGR amplifi- 
cation was carried out in a final volume of 15 \A by reaction 
using equal amounts of cDNAs (2 [A - 200 ng/[A) as tem- 
plate, 0,2 (iM of each primer and 7,5 \i\ of Maxima SYBR 
Green/ROX qPGR master Mix (fermentas, USA) at the 
following conditions: 50°G for 2 min, 95°G for 10 min, 
45 cycles of 95°C for 2 min, 62°G for 30 seg, 72°G for 
30 seg. Data was collected during extension fase. Three 
independent qPGR reactions were performed for final 
quantification. 



Expression levels of GAPDH were used as endogenous 
control. Relative gene expression was calculated using 
the 2~^^^^ method (where GT is threshold cycle) [50]. 
The Pearson correlation coefficient of linear regression 
from 18 pairs of microarray/qPGR expression ratios was 
calculated to validate the qPGR analysis. 
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